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We present new ages and abundance measurements for the pre-main sequence 
star PZ Tel. PZ Tel was recently found to host a young and low-mass companion. 
Such companions, whether they are brown dwarf or planetary, can attain benchmark 
O \ status by detailed study of the properties of the primary, and then evolutionary and 

bulk characteristics can be inferred for the companion. Using FEROS spectra we have 
^ ■ measured atomic abundances (e.g. Fe and Li) and chromospheric activity for PZ Tel 

" and used these to obtain metallicity and age estimates for the companion. We have 

also determined the age independently using the latest evolutionary models. We find 
CSJ ■ PZ Tel to be a rapidly rotating (wsini =73±5kms~^), approximately solar metallicity 

> I star (logA^(Fe)=-4.37 dex or [Fe/H]=0.05 dex). We measure a NLTE lithium abun- 

dance of logAr(Li)=3.1±0.1 dex, which from depletion models gives rise to an age 
of Myrs for the system. Our measured chromospheric activity (logi?^j^ of -4.12) 
returns an age of 26±2 Myrs, as does fitting pre-main sequence evolutionary tracks 
(Te«oi=22±3 Myrs), both of these are in disagreement with the lithium age. We spec- 
ulate on reasons for this difference and introduce new models for lithium depletion 
that incorporate both rotation and magnetic field affects. We also synthesize solar, 
metal-poor and metal-rich substellar evolutionary models to better determine the 
bulk properties of PZ Tel B, showing that PZ Tel B is probably more massive than 
previous estimates, meaning the companion is not a giant exoplanet, even though 
a planetary-like formation origin can go some way to describing the distribution of 
^ . . benchmark binaries currently known. We show how PZ Tel B compares to other cur- 

■ rently known age and metallicity benchmark systems and try to empirically test the 

effects of dust opacity as a function of metallicity on the near infrared colours of brown 
dwarfs. Current models suggest that in the near infrared observations are more sensi- 
tive to low-mass companions orbiting more metal-rich stars. We also look for trends 
between infrared photometry and metallicity amongst a growing population of sub- 
stellar benchmark objects, and identify the need for more data in mass-age-metallicity 
parameter space. 
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1 INTRODUCTION 

PZ Telescopii (HD174429), more commonly known as PZ 
Tel, is classed as a G9IV star and is found to be a vari- 
able star of the BY Draconis type. The star has been stud- 
* E-mail: jjenkinsQdas.uchile.cl ied by many authors since it is a nearby and bright pre- 
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main sequence star. In particular, PZ Tel has been a major 
focus for wor k relating to stellar activi t y and diffe r ential 
rotation (e.g. [Bodp fc Hearnshawl 1 19831 : llnnis et~aL 1 1 19881 : 
iBarnes et al.ll200C ). it is thought to be part of the P Pictori s 
Moving Group, hereafter BPMG, l|Zuckerman et al.|[200ll ). 
and it has been shown to host a debris disk of remnant for- 
mation material of around 0.3 Lunar-masses that spans a 
radius of 165 AU, wit h the inner edge l ocated only 35 AU 
from the central star l|Rebull et al.|[200^ . 

More recently, PZ Tel has been found to be part of 
a binary system w ith a low-mass brown dwarf companion 
l|Biller et al.ll2010t ). Biller et al. used the Gemini-NICI adap- 
tive optics system to discover the faint and co-moving com- 
panion PZ Tel B and from their analysis of its colours and 
luminosity, they found it to have a mass of 36±6Mj, one 
of the lowest mass binary companio n s yet detected. Soon 
after this discovery, iMugrauer et al.l (|20ld ) confirmed the 
existence of PZ Tel B using observations made with the 
VLT-NACO system. They find a mass for the companion of 
28l4^Mj, assuming an age of 12l4Myrs and using the evo- 
lutionary mo dels of lChabrier et alll|2000l ') and lBaraffe et all 
l|2002l , I2OO3I ). in good agreement with the mass found by 
Biller et al. The physical separation of the pair is found to 
be ~16 AU at present, locating the companion within the 
inner edge of the debris disk, which was probably cleared 
out by the formation and evolution of PZ Tel B. 

Binary systems that contain a brown dwarf companion 
and a host star that can give rise to robust evolutionary or 
physical parameters are very useful tools as calibrators for 
modeling the physics of substellar atmospheres. Main se- 
quence stars as hosts can provide very precise m etallicities 
for bro wn dwarf models to benchmark their efforts. [Liu et al.l 
(|2007l ) imaged a faint T-dwarf companion to the exoplanet 
host star HD3651A, allowing a precise metallicity to be as- 
sumed for the T-dwarf. However, robust age estimates for 
field main sequence stars are generally difficult to acquire 
since the evolutionary tracks tend to converge on the main 
sequence. Additionally, main sequence field stars tend to be 
very old, which makes it difficult to ima ge low-mass compan - 
ions that are even known to exist (see lJenkins et ahlboicf ). 
However, a number of surveys are und erway to detect mor e 
brown dwarf binary companions (sec Pinficld ct al. '2005*), 
which have led to recent discoveries (Day- Jones ct al. 2008 : 
IZhang et al.l I2OI0I ) including the discovery of the first T- 
dwarf companion orbi ting any white dwarf primary star 
jPav- Jones et al.ll201ll '). 

1.1 Planetary Origin? 

Another possibility is that PZ Tel B is actually a directly 
imaged exoplanet. 3 0% of disk material c an be locked up 
into forming planets (|Mordasini et al.1l2008l 'l. The maximum 
stable disk mass is given by Md = 0.1 x M©. Therefore, 
the maximum mass of material that can be processed into 
forming rocky cores is given by Mz = Percentage of material 
X Md X Z X (3.3x10^). For instance, a solar mass star can 
contain Mz = 0.3 x 0.1 x 1.0 x 0.02 x 3.3x10^^ ~ 200Mg) 
of rocky material that can form planetesimals. However, for 
PZ Tel A, the disk mass Md would have been 0.10 x 1.16 
~ O.I2M0, and our metallicity measurement gives rise to a 
solar Z fraction of ~2%. Inputing these numbers means that 
the remnant disk around PZ Tel A would have contained 



0.3 X 1.2 X 0.10 X 0.02 x 3.3x10^^ ~ 240Ma) of material 
that could be processed into forming cores. Given that the 
current PZ Tel debris disk only contains 0.3 Lunar-masses of 
material, the majority of the mass of metals could have gone 
into forming a huge core of a few hundred Earth-masses, 
large enough to form a planet like the observed companion 
PZ Tel B; or in this case PZ Tel b. 

This scenario is in agreement with t he latest population 
synthesis models (|Mordasini et al.1l2009l ) where very massive 
planets (20-40Mj or so) can f orm through cor e accre tion, 
but only in a handful of cases. ISahlmann eTall l|2010l ) also 
suggest the dividing line between massive planetary com- 
panions and brown dwarf companions resides in the range 
between 25—45 Mj's. The current mass estimates for the 
PZ Tel companion would place it well within this region. It 
stands to reason therefore that the PZ Tel system could be 
the first system where we have directly imaged one of these 
new extreme- Jovian planets. 



2 OBSERVATIONS AND METHODOLOGY 

2.1 FEROS Spectra 

Two PZ Tel spectra were observed using the Fibre-fed Ex- 
tended Range Optical Spectrograph fFEROS: 'Kauf er et al.l 
|199a) on the MPG/ESO - 2.2m telescope in 2007 and 2010. 
S/N ratios of over 100 in the continuum at 7500A and 
~50-60 at the Ca 11 HK lines (3955A) were obtained and 
FEROS maintains a resolving power of R ~48'000. The 
reduction procedure fo r all spectra is described in more 
detail in Ijenkins et al] (|2008l . I2OIII ). All data were debi- 
ased, fiatfielded, had the scattered-light removed, optimally 
extracted, and had the blaze function removed using the 
FEROS pipeline and a number of Starlink procedures. 

2.2 Metallicity Determination 

Our metallicities were determi ned using a simi l ar me thod 
to that explained in detail in IPavlenko et al.l (|201ll ) ex- 
cept with a few modifications to take into consideration 
the high rotational velocity {vsmi ) of PZ Tel, a method 
that we will outline in a companion paper jlvanvuk et al.l 
I2OIII). We first st a rted with a high resolution and high 
S/N iKurucz et al] l| 19841 ) solar spectrum as our template 
star and then broadened this spectrum to 70 kms~^ to 
closely match t he spectrum of PZ Tel. We used WITA6 
(|Pavlenkd 1 19971 ) to synthesize all spectra by computing 
plane-parall el, and self con sistent, model atmospheres us- 
ing SAM12 (|Pavlenkdl2003l '). along wit h a list of atomic and 
molecular lines dra wn from VALD2 ( Kupka et al.1 ,1999) with 
some updates from lYakovina et alj (j201ll 'l. We then deter- 
mined the best lines/regions in the broadened solar spec- 
trum that would give rise to the well known solar metallic- 
ity value using our spectral synthesis fitting procedure. Once 
these lines were selected, we then select the same regions in 
the PZ Tel spectrum and perform the same analysis as we 
did on the solar spectrum to get the iron abundance of PZ 
Tel. 

Fig. [T] shows the results from our fitting procedure. We 
show best fits in one of our regions of interest, both over 
a wide spectral range (top panel) and a zoomed in region 
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Figure 1. A subsection of tiic FEROS observed spectrum of PZ 
Tel (red) is sliown in a part of our region of interest. The top 
panel shows a wide region and highlights a number of iron lines, 
the lithium lines and other species like aluminium, calcium, and 
silicon. Overplotted in green is our best fit synthetic spectrum 
with the input parameters shown in the key and the blue model 
is a comparison model of higher abundance. The lower panel is a 
zoomed in plot of the lithium lines and shows the same two LTE 
model spectra as in the top panel. 



around the lithium lines (lower panel). In these plots we see 
two LTE model spectra fits (green and blue), with difTer- 
ent abundances, against the observed spectrum of PZ Tel 
(red). Also shown is the model atmosphere we have used 
to synthesize these spectra. In this region we also take into 
account lines of CN and aluminium. It can be seen that the 
models fit the data well, considering the difficulties that one 
encounters by processes such as continuum fitting rotation- 
ally broadened spectra like these, and so we believe that our 
measured metallicity is robust. 

Interestingly, we find the Fe abundance of PZ Tel to be 
log7V(Fe)=-4.37±0.06 dex ([Fe/H]=+0.05±0.20 dex) mean- 
ing it is a solar/slightly metal-rich, young star, a result that 
is not cons istent with most previous metallicity measure- 
ments (e.g. iRocha-Pinto fc M acicl 1998). However, we have 
tested our method on a model young Sun, fit to more than 
one Fe line, and used a methodology we have tested on 
Sun-like stars in the past that was shown to be insensitive 
to sm all variations in the input parameters like T^s and 
log g (IPavlenko et al.ll201ll ). Therefore, our value should be 
more robust than those measured previously. 



Table 1. Calculated metallicity values for differen t temperatures 
around the literature value of 5238 K found by iRandich et al.l 
lll993l) . 



Tea [K] 


logAf(Fc) [dex] 


[Fe/H] [dex] 


4938 


-4.60±0.07 


-0.18 


5038 


-4.53±0.07 


-0.11 


5138 


-4.49±0.07 


-0.07 


5238 


-4.43±0.06 


-0.01 


5338 


-4.37±0.06 


-1-0.05 


5438 


-4.30±0.05 


-1-0.12 


5538 


-4.25±0.05 


-1-0.17 



Table [T] shows the sensitivity of our results over a range 
of Tcff values running from ±300K around the measure d 
Toff of PZ Tel A in the literature (e.g. iRandich et al.l[l993 ). 
The abundances range from -4.60 up to -4.25 in log A''(Fe), 
which relates to a range in [Fe/H] of -0.18 up to -1-0.17 dex. 
This analysis shows that over a span of 600K, we find 
changes in the metallicity of ±0.18 dex, standard deviation 
of ±0.13, and we also see from the table that the uncer- 
tainties on the individual fits are similar across all TcS 's, 
with a little lower uncertainties for the higher T^s models. 
We studied the relationship between the excitation energy 
of the individual lines against their measured abundance, 
and a trend may be present in the data at a Toff of 5238K, 
however we note that most of our lines have similar excita- 
tion energies and only one line had a significantly lower one, 
which determines the slope of the fit. Yet the slope did ap- 
pear to correlate with changing Tes as expected. Therefore, 
we believe the Teff of PZ Tel A to be slightly hotter than the 
previous estimates of 5238K, and so we adopt 5338±200K, 
giving rise to our measured abundance of ±0.05±0.2 dex, in 
good agreement with a solar metallicity for the star. A hot- 
ter Toff can also found from the evolutionar y track s if the 
conversion tabl e from Kenvon fc HartmannI (jl995l ) is em- 
ployed with the lSiess et al. (2000) evolutionary models (see 
next section). 

Finally, given that our Fe lines are located in the opti- 
cal part of the spectrum, and PZ Tel B is 5.04 magnitudes 
fainter than PZ Tel A in the Ks-band, where low-mass brown 
dwarfs emit much more flux than in the optical, we can be 
fairly sure that changes in the metallicity measurements for 
PZ Tel A between different authors is not due to varying 
levels of contamination from the secondary component. 



3 AGE ESTIMATIONS 
3.1 Evolutionary Age 

We first derive the age of PZ Tel A by plotting the star 
on a n HR-diagram us i ng th e late st Hipp arcos data available 
from lPerrvman et al.l ([1993) and Ivan Le cuwcn ( 2007 ). From 
Hipparcos we find the star to have a. B — V colour index of 
0.784 and, with a parallax of 19.42±0.98mas, we obtain a 
distance of 51.49±2.60pc. 

We find that PZ Tel A is most likely a pre-main se- 
quence star, yet to reach the zero age mai n sequence. We 
interpolate its position onto the CESAM (|Marques et al.l 
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Table 2. PZ Tel calculated values 



Star V B-V T^a log g 

(mags) (K) (dex) Mq Kq 

HIP92680 A 8.43 0.784 5338±200 4.41±0.10 1.13±0.03 1.23±0.04 



[Fc/H] log-^HK '°g ^(Li) Te„o; Tgyro ri,i 

(dex) (dex) (kms"^) dex (Myrs) (Myrs) (Myrs) 



+0.05±0.20 -4.12±0.06 72±5 3.1±0.1 22.34±3.13 26.23±1.62 T 



l2008h and lSiess et aD l)200d ) isochrones and isomass tracks 
in a similar manner to that in IJenkins et all l|2009l ). The CE- 
SAM models give rise to a mass of 1.14±O.O3M0 , a surface 
gravity (log g ) of 4.41±0.10 dex, radius of 1.26±0.04Rq and 
an age of 22.17±2.05Myrs. These are in good agreement with 
the values from the Siess et al. evolutionary models which 
also take into consideration the metallicity of the star. The 
mass, radius and age from these models are l.llzbO.OSM© , 
1.20±0.03Rq and 22.51±2.36Myrs, respectively. The final 
values are listed in Table [2 



3.2 Gyrochronological Age 

We extract the FEROS chromospheric logi?HK activity at 
two epochs for PZ Tel, once in 2007 and again in 2010. This 
gives a fairly robust activity measurement however the ac- 
tivity cycle is expected to be fairly large therefore we may 
still be miss representing the true mean activity level. We 
measured these activities following the same procedure ex- 
plained in lJenkins et al.l (|2006l . 120081 '). Briefiy, the Ca II HK 
lines, located at 3968. 470A and 3933. 668A respectively, were 
filtered through two triangular bandpasses centered on their 
line cores and with FWHM's of 1.09A. We then compare the 
ratio of these filtered fiuxes to another two filtered square 
bandpass regions in the surrounding continuum, centered at 
389lA and 400lA, labelled the V and R bands respectively. 
This ratio is highlighted in Eqn.[l] where the Ni is the inte- 
grated fiuxes in each filtered bandpass region. 



Nferos ~ 



N„ + Nk 
Nv + Nr 



(1) 



This ratio is then normalised to the bolometric lumi- 
nosity of the host star to extract the chromosp heric part of 
the sp ectral light and using the relations from lNoves et all 
l|l984h we arrive at the final logJ^JjK activity index. For PZ 
Tel we find logiinK activities of -4.16 and -4.07 dex for the 
2007 and 2010 data respectively. The mean of these mea- 
surements (-4.12zb0.06 dex) is used to derive its age. The 
difference between the two values highlights the young and 
active nature of the star, given the uncertainties are only at 
the level of ±0.02 dex for these observations. 

In order to obtain an age estimate from the activ- 
ity of PZ Tel we use the latest age-activity relations from 
iMamaiek fc Hillenbranll (|2008l ). The Mamajek & Hillen- 
brand relation gives rise to a mean age of 26.23±1.62 Myrs 
for PZ Tel, with a range between 17—40 Myrs for the two 
individual measurements. The uncertainties on the age es- 
timation were taken using the published scatter around the 
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Figure 2. Top panel shows the differences in measured equiv- 
alent width of lithium by adopting LTE (solid curves) or NLTE 
(dashed curves) . The horizontal arrow shows the measured values 
for PZ Tel A for two plausible T^ff values. The lower panel shows 
lithium depletion as a function of age taken from the ISiess et al.l 
pre- main sequence models for a 1.15M0 star. The cross 
represents our measured value for PZ Tel A and the key in top 
right highlights the measured age, lithium abundance, and metal 
fraction employed. 



age-activity fit, along with the uncertainties on the mean 
age derived from the two observations. This age estimation 
is in very good agreement with the age we have derived from 
the isochrone fitting procedure above. 
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3.3 Lithium Age 

We also measure an age for PZ Tel through model fitting the 
Li lines at 6708 A (as shown in Fig. [T|. Strong lithium ab- 
sorption is another indicator of youth in stars since lithium 
is destroyed at temperatures above around 2.5xiO®K, and 
therefore this element is depleted through time in the in- 
teriors of stars. The depletion of lithium is dependent on a 
number of factors that influence the convective envelopes of 
stars and the mixi ng processes and convection therein (see 
iPinsonneauldl 19971 ). Yet lithium can still be used as a good 
indicator of the age of young stars, even though some stud- 
ies have shown that evolutionary models may under-predict 
the lithium depletion, giving rise to systematically ol der ages 
l|Zuckerman et al.ll200ll : IWhite fc HiUenbrandllioOsI ). 

In the lower panel of Fig. [1] we see our best flt syn- 
thetic spectra to the lithium region. The models we show 
fit the observed spectrum well and we find the best fit to 
be for a wsini of 73 kms~^ and a LTE lithium abundance 
of log 7VLTE(Li) 3.3. We then follow the NLTE co r nputa - 
tional procedure explained in iPavlenko fc Magazzul (|l996l ) 
to arrive at our final value for the lithium abundance of PZ 
Tel of log A^NLTE(Li) 3.1±0.1 dex, in good agreement with 
that found by previou s authors (e.g. iRandich et al.l 1 19931 : 
ISoderblom et al.il 19981 ). 

We highlight that our value is drawn from a NLTE anal- 
ysis, since the metallicity we measure for PZ Tel A was found 
under the assumption of LTE. The top panel in Fig.[2]shows 
the difference in equivalent width of the lithium line as a 
function of measured abundance. The LTE curves have a 
shallower gradient at a higher log A''(Li) than the NLTE 
curves and both cross at a log A'^(Li) of around 2.6-2.8. We 
show two different temperature tracks and we see that in 
general, models of higher TbS provide higher lithium abun- 
dances. Our LTE spectral synthesis provides the best fit Li 
abundance of 3.3 for the lithium abundance, shown by the 
horizontal arrow in the figure, but this value drops to our 
measured value of 3.1 when examined in NLTE. 

Fig. [5] (lowe r panel) shows the lithium depletion evolu- 
tion models from lSiess et al.l (|2000l ) for the measured mass of 
PZ Tel A (I.I5M0). Lithium models generally are for non- 
rotating stars, whereas young and fast rotators like PZ Tel 
induce processes difficult to model accurately. Our measured 
lithium abundance for PZ Tel is highlighted on this plot by 
the cross. We also show the age and abundance value in the 
key at top right of the plot. We measure a lithium depletion 
age for PZ Tel of only 7I2 Myrs, younger than both the 
chromospheric age and the evolutionary model age. Also, 
the chromospheric age and evolutionary model fitting ages 
are far older than the depletion timescale for a star at this 
Tcff by a factor of two or so. Such discrepancies require expla- 
nation, particularly since the evolution of Teff as a function 
of age from the same models yield an age of 23 Myrs. 



3.3.1 Rotation, Magnetic Fields and Accretion 

As mentioned, the depletion models are non-rotating mod- 
els, whereas PZ Tel A is a fast rotating star. We could envis- 
age th at significan t rotation could drive powerful magnetic 
fields (|Kraftlll967l ). which could inhibit the convective mo- 
tions within the young stellar atmosphere. Inhibiting this 
motion could lead to reduced lithium depletion in the inte- 
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Figure 3. Three 1.15M0 models of litliium depletion for a stan- 
dard, non-rotating, and non-magnetic star (solid), a magnetic star 
(dashed), and a fast rotating star (dotted). 



rior of the star (jMartin et al.|[l993 ). In addition, these mag- 
netic fields could also reduce the level of differential rota- 
tion in these stars, yieldin g less shearing and again further 
reduced lithium depletion. iBaxnes et al.l (|2000l ) has studied 
the level of differential rotation in PZ Tel A and found the 
surface shear to be similar to the solar shear, possibly sup- 
porting this scenario, but this is also simi lar to other young 
stars like Speedy Mic (|Barnes et al.ll200ll ). 

We have investigated this issue using new models cur- 
rently unde r construction based around the explanations 
presented in lChabrier et akj l|2007h who introduced activity 
effects in low-mass stellar models by considering two param- 
eters: the spot blocking factor (/3) and the modification of 
the mixing length parameter (a). This study shows that the 
effects of the a and /? parameters are degenerate, i.e. the 
properties of any given system can be reproduced by modi- 
fying any of the two or both. 

The authors also showed that the introduction of 
rotation and / or m agnetic field effects on the models of 
iBaraffe et al.l l|l998l ) could explain the discrepancies between 
the observed and theoretically predicted mass-radius rela- 
tionship of eclipsing binaries. Specifically, they presented 
two scenarios considering: (1) that the effect of magnetic 
fields and rotation alter the efficiency of convective energy 
transport, which can be modelled by setting the mixing 
length parameter (Mixing Length Theory) to lower values 
than those used for solar models; and (2) that the stellar 
magnetic activity present on these objects can be associated 
with the appearance of dark spots covering the radiative sur- 
face, modeled by the P factor indicating the percentage of 
the stellar surface c overed by spots . 

The results of IChabrier et all l|2007l ) show that both 
these scenarios predict larger radii than standard stellar 
models, but while the effects of spots are significant over 
the entire low-mass domain, the effect on convection is rel- 
atively small for fully convective stars. Moreover, modified 
evolutionary models present cooler central temperatures, af- 
fecting burning rates of light elements. Specifically, lithium 
is depleted more slowly in models where the convection is in- 
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hibited or where spots are present. If the lithium is burning 
at different rates, the age derivation for stellar or substellar 
objects needs to be revised due to the effects of activity or 
rotation in t his calculation. Work is under way to explore 
these effects (|Gallardo et al.ll2011^ . 

Fig. [3] shows preliminary results for the Li evolution 
of a I.I5M0 object for a standard model (solid line, spot- 
free and Qf=1.9), and two "active" models: short-dash line 
(a=1.9 and /3=0.5 i.e., 50 % surface coverage by spots) and 
dotted line (fast rotation, a=0.5 and spot-free). As we can 
see from the plot, the "active" models significantly decrease 
the rate of lithium depletion, through inhibiting the convec- 
tive motions in the stellar interior. Thus, rotation and/or 
magnetic fields (through spot coverage) could be present on 
our target and could give rise to the high Li abundance we 
measure for PZ Tel A, meaning a younger age from a non- 
rotating model could be expected. 

An additional possibility is that the lithium depletion is 
not being affected at all by magnetic fields and high rotation, 
but that the level of lithium is being replenished somehow 
in the atmosphere of PZ Tel A. Given that we find PZ Tel A 
to have around solar metallicity we expect there was a lot 
of dust in the proto-planetary disk and hence planetesimal 
formation was a strong possibility. Accretion of planetesi- 
mals could significantly replenish the level of lithium in the 
atmosphere of PZ Tel A, giving rise to the measured value 
of 3 .1 that we currently find fo r this star. 

Ilsraelian et al.l(|200ll . l2003l ) suggest that the planet host 
star HD82943 has engulfed a planet of around 2Mj , leading 
to an enhancement of the ®Li isotope up to the measured 
value of 4.5x10** nuclei, therefore a similar process could 
explain the overabundance of lithium in PZ Tel A's atmo- 
sphere. For such young and rapidly rotating stars it is diffi- 
cult to conclude this is the case since the uncertainty on the 
lithium abunda nce is large. 

In addition, iBaraffe fc ChabrieJ (|2010l ) have shown that 
episodic accretion onto the star can significantly affect the 
lithium abundance, adding another source of uncertainty. 
Given these results and conclusions, we do not use the age 
derived from the lithium abundance in our final mean age 
for the PZ Tel system. 

Taken all together, we get an average age of 24±3 Myrs 
for PZ Tel, which is slightly higher than the previous esti- 
mates g iven that PZ Tel is tho ught to be a member of the 
BPMG IZuckerman et al.ll200il '). but only at the level of la. 
However, this could indicate that PZ Tel is not a bonafide 
member of the BPMG, or that the true age of the BPMG 
is actually significantly older than the I2I4 Myr age esti- 
mated by Zuckerman et al. It is important to remember 
though that the uncertainties quoted on the age here are 
formal and are almost certainly a lower limit, with the true 
uncertainty larger than this quoted value. 



4 DISCUSSION 

4.1 Moving Group Member? 

The age and metallicity for PZ Tel A measured in this work 
is slightly different from those previously published in the 
literature for this star. The age we measure is greater than 
the age derived by other authors of I2I4 Myrs, based pri- 
marily on the notion that the PZ Tel system is part of the 



BPMG, only at the Icr level. We find a mean age for PZ 
Tel A of ~24±3 Myrs, based on three different methods, 
two of which are in good agreement. It is difficult to say 
whether this age can be transferred to the entire BPMG, 
if we assume that the PZ Tel system is a bonafide member 
of this kin ematic group. However, given the new Hipparcos 
reduction (|van Leeuwenli2007 ) has altered the evolutionary 
status of PZ Tel A and the star's high rotational velocity 
(which makes the measurement of accurate radial velocities 
very challenging) , it may be necessary to reassess both if the 
PZ Tel system is an actual member of the BPMG and also 
the age of the BPMG itself. On the other hand, this result is 
more in agreement with the age es timated for the BPMG by 
iBarrado v Navascues et al.l (|l999h of 20±10 Myrs, based on 
their discovery of three new M-dwarf members to this mov- 
ing group, indicating that the BPMG is indeed older than 
12Myrs. 

The metallicity ([Fe/H]) we measure for PZ Tel A is 
found to be -|-0.05±0.20 dex, more than 0. 30 dex higher 
than the earliest estimates for this star (|Randich et al.l 
ll993l : lRocha-Pinto fc Maciellll998l ). how ever more in agree- 
ment with the metallicity quoted later bv lRocha-Pinto et al.l 
(l2000l ) of 0.16 dex. It is interesting to note that this later 
measurement by Rocha-Pinto et al. was made photometri- 
cally, where as the most recent pho tometric metallicity es- 
timate from iHolmberg et al.l l|2007l ) provides a very metal- 
poor photospheric metallicity abundance of -0.50 dex. This 
highlights the uncertainties inherent when using photo- 
spheric colours and magnitudes to estimate the metallic- 
ity of stellar atmospheres, particularly for young stars that 
show evidence for variability. We argue that our metallicity 
should be more robust than past measurements, given the 
steps we explained above to ensure we selected the best Fe 
lines in our analysis. Indeed, lower metallicity values were 
found when we randomly selected all lines, but this sys- 
tematically affected the solar metallicity also by shifting it 
to more metal-poor values. This may be the reason for the 
metal-poor nature found by Randich et al., along with their 
use of older model atmospheres and more incomplete spec- 
tral line lists. 

The solar/metal-rich nature of the PZ Tel system agrees 
well with the notion that young stars are formed in more 
metal-rich environments given enrichment of the interstel- 
lar medium by past supernovae explosions. In addition, it 
also ties in with the finding of a debris disk around PZ Tel, 
which would indicate there was an abundance of metals in 
the remnant disk of material leftover by the formation of 
PZ Tel A, even though the current disk is estimated to con - 
tain only 0.3 Lunar- masses of material (|Rebull et al.llioO^ '). 
We searched the literature for spectroscopic metallicity val- 
ues measured for other BPMG members, yet none wer e 
found, pa rticularly in the works of Zuck erman et al. (200l|), 



iFeigelson et al.1 (|2006l ) and lOrtega eta l. (2009). Therefore, 
if PZ Tel A is a BPMG member then this indicates that the 
BPMG is a cluster of around solar metallicity, if all members 
do indeed share simil ar metallicities. Wor k is underway to 
explore this question l|lvanvuk et al.|[2011^ . 

4.2 Evolutionary Model Testing 

Fig. |4] (upper panel) shows the position of PZ Tel B 
on a colour-magnitude diagram. The evolutionary calcu- 
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lations for such low mass objects are based on the Lyon 
stellar evolution c o de w ith in put physics d escrib ed in 
IChabrier fc Baraffd (| 19971 ) and iBaraffe et all l| 19981 ). The 
outer boundary conditions between the interior and the 
non-grey atmosphe r e pro files are presented in detail in 
IChabrier fc Baraffd (120001') as well as the treatment of dust 
in the atmosphere in Allard et al.l (|200ll ). We show them 
overplotted for metallicities of 0.00, +0.05 (metallicity of PZ 
Tel A), and +0.50 dex, ages between 1-50 Myrs, and masses 
from 10-lOOMj. The position of PZ Tel B is represented by 
the open circle with its associated uncertainties, which in 
colour are fairly su bstantial. The valu e is a ctually the mean 
values taken from iBiller et al.l (|201(J ) and iMugrauer et al.l 
l|201(tl . 

Firstly, we find that for young ages and such low-masses, 
increasing the metallicity from solar to +0.5 dex causes an 
increase in effective temperature at the level of ~250K, at 
the upper mass limits for brown dwarfs (~70Mj). Around 
PZ Tel B we find this value is lower, at the level of ~150K. 
This highlights that increasing the metallicity of low-mass, 
young substellar objects has a significant effect on the bulk 
thermal properties of the system, particularly for a signifi- 
cant increase in metallicity at the +0.5 dex level. 

Given the rather large uncertainties in the J — H colour 
for PZ Tel B, it does not facilitate a useful diagnostic at 
present to distinguish between different evolutionary tracks. 
To la PZ Tel B is in agreement with all the model tracks 
here. Even so, currently the evolutionary tracks are not ro- 
bustly tied to observations of substellar objects with well 
determined metallicities and gravities and therefore many 
more benchmark objects are needed to facilitate model con- 
straints. 

The lower panel in Fig. [3] shows the position of PZ 
Tel B as a function of our newly dete rmined age, against 
the mean Mj mag n itude from both iBiller et al.l (|2010l ) 
and lMugrauer et al. (|2010l ) (filled circle). The evolutionary 
tracks are again plotted, with the solid curves representing 
the metal-rich model of +0.05 dex and the dashed curves 
representing the solar metallicity models. We can firstly see 
that there is not a significant difference between the two 
evolutionary tracks, however for rigour we have used the 
+0.05 dex models in all our calculations. We find PZ Tel B 
to have a mass of 62+2 Mj. This value is around twice the 
quoted absolute mass values given in Biller et al. and Mu- 
grauer et al., however we note again that the uncertainties 
on the age are formal and are certainly a lower limit and 
we quote a more realistic uncertainty from our Monte Carlo 
simulations that we explain below where we consider not 
just one fixed point in the parameter space but analyse our 
overall range of values. The TcS and log g measurements we 
obtain are 2987+100 K and 4.78+0.10 dex, respectively. 

In the figure we also show the position of PZ Tel B given 
the Zuckerman et al. estimates for the age to the BPMG 
used in both Biller et al. and Mugrauer et al. to obtain their 
mass estimates for PZ Tel B (open square). The mass we 
derive from the models using this age estimate is 32lg'*Mj, 
which runs from ~24-56Mj. Therefore, although the mass 
estimate we obtain given our new age measurements for the 
PZ Tel system is almost twice this value, the significance 
of the difference is only at the level of around 1.2a. The 
Teff and log g measurements we obtain are 2377I34QK and 
4.68lQ'j7dex, respectively. 
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Figure 4. The top panel shows a set of low-mass substellar evolu- 
tionary tracks on a colour-magnitude diagram. The key highlights 
the ages as a function of colour and also we show solar metallic- 
ity models (dashed curves), our measured PZ Tel A metallicity 
of +0.05 dex models (solid curves), and super-solar metallicity 
(+0.50 dex) models (dot-dashed curves). PZ Tel B is represented 
by the open circle. The lower panel shows the change in abso- 
lute J-band magnitude as a function of age, for different substel- 
lar masses. The solid and dashed curves represent the measured 
metallicity and solar metallicity models, respectively. The filled 
circle is the position of PZ Tel B, given our age estimate for the 
syste m, whereas the op e n squ are represents the age for the BPMG 
from IZuckerman et al.l 1 I2OOII) . assuming the PZ Tel system is a 
true member of the group. 



We note that ISimon fc Schaeferl l|201ll ') have recently 
measured the diameters of two BPMG members HIP560 and 
HIP21547 using interferometric methods, and they find that 
evolutionary models would systematically overestimate the 
masses of these companions by around 0.2 M0 and hence 
the ages are older by around 5 Myrs, both due to the effects 
of gravitational darkening. They find an age of 13+2 Myrs 
for the BPMG. If we take these systematic offsets and ap- 
ply them to PZ Tel A, we get a mass and age for this star 
of 1.03+0.04 Mq and 19.83+2.94 Myrs. This gives rise to 
a mass, %b , and log g of 57tlo ^J' 2923+150 K, and 
4.78+0.10 dex for PZ Tel B, stiU placing it securely above 
the planetary limit. 



4.3 Monte Carlo Analysis 

To better investigate the distribution of possible bulk prop- 
erties for PZ Tel B given our measured input parameters for 
PZ Tel A, we perform a Monte Carlo (MC) analysis on the 
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Figure 5. The best fit parameters for mass (top), T^ff (middle), 
and log g (bottom) as a function of both metallicity (left column) 
and age (right column) for each of our 10000 random realisations 
of the possible values of metallicity, age, and J-band magnitude 
for PZ Tel B. 



data. To do this we randomly varied the input parameters 
age, metaUicity, and J-band magnitude by their associated 
uncertainties, assuming a Gaussian model. For each random 
realisation we then reinterpolate the evolutionary models 
and determine the best fit mass, T^b , and log g . 

We performed 10000 random realisations in our simu- 
lation and the distributions for each are shown in Fig. [S] 
In mass we see a well defined population as a function of 
metallicity, with an indication of a slight trend whereby 
metal-richness leads to less massive companions. We also 
ran the same tests in the H and Ks bands to test if this 
trend towards lower-mass companions in metal-rich systems 
was found across the three near infrared bands, and not 
only the J-band. We found the trend was apparent in all 
three bands, therefore current models suggest that in the 
near infrared, one is more sensitive to lower-mass compan- 
ions around more metal-rich stars. As a function of age we 
see the well understood trend between age and mass. Taken 
together, higher metallicity and younger systems will pro- 
duce lower-mass substellar objects, and hence future sur- 
veys could target the most metal-rich young stars to try 
to bias their searches towards lower-mass objects. This sce- 
nario actually agrees with the prevalence of gas gian t plan - 
ets towards more metal-rich stars ([Fischer fc Valentill2005l ), 
particularly when compared to field binary stars. 

The distribution of both Tcff and log g show a num- 
ber of features of this simulation, such as the truncation of 
the model boundaries, along with very non-symmetrical dis- 
tributions. The Tefi distributions generally follow the same 
trends as the mass distributions for both metallicity and age, 
since mass and Teff are heavily correlated at young ages for 
low-mass substellar brown dwarfs. The distributions have 
parabolic boundaries, centered around the measured metal- 
licity and age of the PZ Tel system, whereby the cooler mod- 
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Figure 6. Gaussian histogram distribution in mass for the Monte 
Carlo analysis performed to determine the bulk properties of PZ 
Tel B. Also overplotted is the best fit gaussian model to the data, 
along with the measured parameters in the top left. 



els are spread more tightly around the measured values in 
comparison to the hotter models, which are generally found 
across a wide range in both metallicity and age. 

As for log g against metallicity and age, we find a dense 
population of values with a high surface gravity clustered 
around 4.9, and then a more collimated clustering towards 
lower gravities. In the metallicity plane we see that for metal- 
poor values the log g can have many solutions between 4.75 
to 4.9, however the metal-rich models produce a cluster- 
ing towards higher gravities. A similar trend is seen in age, 
where young ages produce higher surface gravities, when 
compared to ages higher than the measured age of the PZ 
Tel system. In fact, there are not many solutions in agree- 
ment with the measured log g from the final position in age 
and metallicity. 

The MC simulations were also used to better define our 
measured uncertainties for the parameters of PZ Tel B. An 
example of this is shown in Fig. [6] where we plot the fre- 
quency distribution of masses in histogram format for the 
entire set of 10000 random realisations. The solid curve over- 
plotted is the best fit gaussian model to the data and we 
show the measured parameters from this best fit model. We 
find the mean of the mass distribution to be 62Mj and the 
standard deviation is 9Mj. Clearly this is much larger than 
the formal uncertainties of ~2Mj quoted above for the mea- 
sured solution and therefore we adopt this uncertainty for 
our mass measurements. We performed a similar analysis for 
the Tcff and log g and quote these uncertainties in Table [S] 



4.4 Benchmark Systems 

Fig. [7] shows the position of PZ Tel B (open circles) with 
respect to other potential age-metallicity benchmark binary 
systems currently known (filled circles) that were taken from 
iDav- Jones et al.] (|201ll ). In the top panel we show metallic- 
ity against system age in Gyrs and the lower panel shows 
metallicity against companion mass, in Jupiter-masses. The 
characteristics of these benchmark systems are shown in Ta- 
ble[3l along with our calculated values for PZ Tel B for direct 
comparison. 

The position of PZ Tel B in metallicity and age space re- 
veals just how this object extends the benchmarks into a new 
part of the parameter space. PZ Tel B is far from any other 



© 2002 RAS, MNRAS 000, [T]-?? 



Age and Abundances of PZ Tel 9 



10.00 
1.00 

0.10 

0.01 
80 

^ 60 

\ 40 



20 - 



®- 



-0.6 



-0.4 



-0.2 0.0 
[Fe/H] 



0.2 



0.4 



Figure 7. The current distribution of age-metallicity benchmark 
binaries. The top panel shows metallicity against age, with sym- 
bol sizes scaled by increasing mass of the companions. The dot- 
ted cross-hairs represent the solar values in metallicity and age. 
The lower panel shows metallicity against mass with the symbol 
sizes scaled by increasing age of the systems. In both panels PZ 
Tel B is shown as the open circle with its associated uncertainties, 
the benchmark systems are represented by the filled circles, and 
the ringed filled circles are the directly imaged planets around 
HR8799. Also shown are formation models for core masses of 1% 
(dashed curve) and 5% (solid curve) of the total mass of the com- 
panions. 



benchmark binary object currently known, potentially giv- 
ing a new insight into the physics of solar metallicity, young 
brown dwarfs, once a spectrum of the object is obtained. The 
symbol size representing each object is scaled to the mass of 
the most massive benchmark object, HD89744 B. The mass 
of PZ Tel B is around the median of the other benchmarks, 
which again means PZ Tel B could be an interesting test 
case for studying low-mass atmospheric physics. However, 
only 77 Cancri B has colours, or a spectral type, approaching 
that of PZ Tel B, even though this object is at the oppo- 
site end of the age distribution. We do note though that 77 
Cancri B is likely an unresolved LT binary and therefore the 
properties for this object may not be accurate. 

In this plot we also show the position of the Sun, in 
both metallicity and age, marked by the dotted lines, where 
the measured values are where these lines cross. This al- 
lows us to see how the benchmarks compare to the Sun. All 
benchmarks are younger than the Sun, which represents a 
strong bias towards younger objects, given they are much 
more luminous for a given mass, with AB Pic B and PZ 
Tel B much younger than the Sun and the other benchmark 
binaries. Also we see that most are fairly close to the solar 
metallicity by around ±0.2 dex, including PZ Tel B, except 
the extremely metal-poor brown dwarf binary, AB Pic B. 
Taken all together, it seems that we have covered a large 
fraction of the benchmark metallicity-age space, except for, 

1) all metallicities with ages older than the Sun (^4 Gyrs), 

2) metallicities less than -0.2 with ages above ~0.1 Gyrs, and 



3) young ages (^1 Gyrs) with solar/super-solar metallicity 
(^0.0 dex), which is where PZ Tel B finds itself positioned. 

The lower panel shows the mass range for the popula- 
tion against metallicity and there is a clear clustering of ob- 
jects with high brown dwarf masses and in more solar/metal- 
rich environments. Higher masses would be favoured as there 
is a strong bias towards such objects, given they appear 
brighter on the sky for a given distance than lower mass 
objects and this can go some way to explaining this mass 
clustering. In this case the symbol sizes are scaled by age, 
from the oldest benchmark rj Cancri B. Interestingly, PZ 
Tel B is found clustered around a number of higher mass 
brown dwarfs, both in metallicity and mass space, however 
it is clearly much younger than the others. PZ Tel B actu- 
ally completes a missing age piece of the evolutionary scale 
in this mass and metallicity parameter space, meaning this 
region is reasonably well sampled in comparison to other re- 
gions. Probably the most important point here is that there 
is a lack of low-mass objects across all metallicities, and 
the three objects with masses less than 60Mj, are found to 
be very young, which if PZ Tel B is really younger than our 
measurement of ~24 Myrs, then it will drop into this regime, 
but again would be a very young object. 

We also note that for these binaries there may also be a 
formation effect at play here too when looking at mass, de- 
pending if the companion has formed through direct grav- 
itational collapse of the remnant disk, or if it has formed 
through a core accretion process, more akin to planet for- 
mation. Core accretion models generally do have difficulty 
forming massive objects, however on rare occasions, models 
do predict such large bodies can form through the planet 
formation process. 

We overplot two core accretion based models, using the 
relationship discussed in § 1.1. We have made some broad as- 
sumptions here, for instance, the companions all have a fixed 
envelope to core mass fraction. The two curves shown repre- 
sent 5% (solid) and 1% (dashed) core mass fractions, and we 
have extended the 1% model up to the brown dwarf bound- 
ary. Physically this may be unrealistic as the core accre- 
tion models have prob lems building such large objects (see 
iMordasini et al.ll2009l ). Fragmentation of the remnant disk 
might be a better way to form these high er mass compan- 
ions (see IStamatellos fc WhitworthI 120091 ) . However, these 
basic models do show that when there is a high fraction of 
metals in a proto-planetary disk, and extremely large cores 
can be formed, if we input realistic core-to-envelope ratios 
of only a few percent we can go a long way to explaining 
the distribution of benchmark binaries in metallicity-mass 
space. 

In both the upper and lower plots the filled circle en- 
cased by the rings mark t he position of the d irectly imaged 
planets HR8799 b, c and d l|Marois et al.ll2008l ). If these wide 
orbiting planets were formed by core accretion, along with 
the benchmark brown dwarf binaries, then it is worth com- 
paring the properties of these planets with the brown dwarfs. 
Of course, again we have the biases, particularly in age, 
where these planets were only able to be imaged with current 
technology since they are so young, and hence bright. How- 
ever, in the metallicity-mass parameter space, these planets 
reside between the two formation models, along with a few 
of the benchmark binaries, including PZ Tel B. 
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Figure 8. A colour-colour diagram showing how PZ Tel B 
(open circle) compares against other benchmark binaries cur- 
rently known (filled circles). Also shown are the planets around 
HR8799 (open diamonds) and a list of free-floating planetary mass 
objects (crosses). The filled circles with rings around them show 
objects with masses and ages close to PZ Tel B. The solid and 
dashed curves are DUSTY and COND evolutionary models re- 
spectively. 



4.5 Benchmark Colours 

4.5.1 J-H vs H-K 

Late-M stars are known to be affected by emergent dust in 
their a tmospheres at Teff ' s below the temperature of PZ 
Tel B (| Jones fc Tsuiil Il997l ). They also exhibit features in 
their near infrared spectra such as water bands around 1.4, 
1.8 and beyond 2.4/im, CO bands at around 2.35/im, FeH 
bands at 0.99 and 1.2/im, strong J-band potassium absorp- 
tion. We show the position of PZ Tel B in comparison to 
our other benchmark systems on a colour-colour diagram 
{H — Ks against J — H ) in Fig. [51 Also sho wn for com- 
parison are the DU STY (jChabrier et al.ll2000l ) and COND 
jBaraffe et aLll2003l ) evolutionary tracks, represented by the 
solid and dashed curves, respectively. 

First of all, we find that PZ Tel B is much bluer in 
both colour bands than the other benchmark binaries we 
have discussed, including the binaries with masses and ages 
near that of PZ Tel B. It is also bluer than some young 
planetary-mass M dwarfs from the Orion Cluster that we 
also highlight (; Weights et all |2009| ). at least bluer in the 
H — Ks colour index. This could be an affect related to the 
difference in metallicity between PZ Tel B and the Orion 
Cluster objects since some estimates for the metallicity o f 
the cluster place it below the solar value (see IO'delll200ll ). 
however later estima tes claim a more solar metallicity value 
l|D'Orazi et al.ll2009l '). 

These low-mass and young free-floating planetary-mass 
objects are useful to interpret PZ Tel B, since they are of 
similar age and T^ff , and they help us to realise that PZ 
Tel B could be an M, L, or T-dwarf or even a planet. Again 
this raises questions on formation scenarios for young and 



low-mass objects, but questions over the true physical nature 
of P Z Tel B will require spectroscopic follow-up. 

iLeggett et all l|200ll ) note that when there is no dust 
present in the atmospheres of M dwarfs, the water bands are 
expected to become deeper and exhibit increasingly steep 
wings with decreasing Teff . However, in the presence of dust, 
the atmosphere is heated and the bands become wider and 
shallower. We have shown that PZ Tel B is not a metal- 
poor companion through association to the host star PZ 
Tel A, and hence PZ Tel B can be expected to host a dusty 
atmosphere. In particular, we predict the water bands for 
this companion to be shallow and broad in comparison to 
more metal-poor field M dwarfs of a similar Toff . 

Another opacity source affecting such cool atmospheres 
is tha t of coUisionally induced hydrogen absorption (Linskvl 
Suppressed K-hand flux has been attributed to this 
and has been used in colour selection criteria to select un- 
usual substellar objects (e.g. [Murray et ai]|2011^ . However, 
since PZ Tel B does not appear to be significantly redder 
in the H — Ks colour band compared with the other young 
M dwarfs we show, it appears that this type of absorption 
is not a significant source of opacity in the atmospheres of 
young and low-mass brown dwarfs like PZ Tel B. Given that 
we only have one data point here at present, this may po- 
tentially provide a future avenue of research. 

Although PZ Tel B is bluer in both colours in com- 
parison with the companion planets of HR8799 (open di- 
amonds), the free-fioating planetary-mass objects exhibit 
similar colours to these planets. This shows that the evo- 
lutionary properties of young low-mass brown dwarfs and 
high-mass planets are similar, validating their use as bench- 
marks to better understand the physics of gas giant planets. 

4.5.2 Metallictty vs Near-IR Photometry 

Since metallicity can have a strong impact on the atmo- 
spheric opacities of cool substellar objects, studying how 
colours evolve as a function of metallicity can allow a 
deeper understanding of the physics and interactions ongo- 
ing within these objects. In Fig. [9] we show the distribution 
of J — Ks and H — Ks colours against changing metallicity. 
Given that substellar colours evolve as a function of time 
and as a function of mass, each object's colour was scaled to 
the age and mass of PZ Tel B using the DUSTY evolutionary 
models ijChabrier et al.ll2000D . This scaling allows us to at- 
tempt to isolate the effects of metallicity on the atmospheric 
properties of these companions. 

At first glance there maybe a trend between the 
metallicity and the broadband colours of these bench- 
mark companions, even when including the planets around 
HR8799 (open diamonds). A possible anti-correlation is 
present, whereby an increase in metallicity tends to de- 
crease the colour index, at least in the metal-rich regime. 
The Pearson rank correlation coefficient for the ,JKs data 
is -0.73, suggestive of a strong trend between these two 
parameters. Such a trend would indicate a turnover is 
present, since the met al-poor L sub-dwarfs also exhibit 
bluer JKs colours (e.g. iBurgasser et al.l [20091 : iLodieu et al.l 
I2OIOI ). If this is the case, then it would provide an expla- 
nation for at least so me of the peculi a r blue L-dwarfs dis- 
cussed in length in iKirkpatrick et al.1 l|2010l ). Objects like 
2M171H-4028 B, which have been shown to have solar 
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Figure 9. The distribution of near infrared colours J — Kg (top) 
and H — Ks (bottom) as a function of increasing metallicity for 
PZ Tel B (open circle), the comparison benchmark binaries (filled 
circles), and the planets orbiting HR8799 (open diamonds). All 
object colours were scaled using evolutionary models to the age 
and mass of PZ Tel B for comparison. Uncertainties are included 
where they were published in the literature. 




-0.6 



-0.4 



-0.2 0.0 
[Fe/H] 



0.2 



0.4 



Figure 10. The near infrared absolute magnitudes against metal- 
licitics for PZ Tel B (open circle), the benchmark binaries (filled 
circles), and the planets around HR8799 (open diamonds). From 
top to bottom we show absolute J, H, and Ks magnitudes, re- 
spectively. The absolute magnitudes of the benchmarks have been 
scaled to the age and mass of PZ Tel B for direct comparison. 



metallicity (|Radigan et al.l |2008| ) and blue colours, would 
naturally be explained if there is a negative trend between 
near infrared colours and metallicity in the more metal-rich 
regime. 

If at least part of the correlation discussed is real, it 
could be more important in the bluer bands than in the red- 
der bands, as it affects the JKs colour more than it affects 
the HKs colour. The Pearson rank coefficient for the HKs 
colour is -0.65, which still indicates a trend, but at a reduced 
level of significance. 

In Fig. 10 we show the absolute J, H, and Ks magni- 
tudes against metallicity to test if any gradients are apparent 
in the absolute fluxes that would support the possibility of 
a decreased level of dust or absorption is present. Indeed, 
we do see a possible anti-correlation between the near in- 
frared photometry and increasing metallicity, with the trend 
appearing a little stronger in the J-band than the redder 
bands. 

However, given that the two metal-poor data points are 
both young and of low-mass there could still be an age and 
mass correlation present which has not been fully corrected 
for due to model uncertainties. The models have known diffi- 
culties in reproducing the bulk properties of young and low- 
mass brown dwarfs. To test this further it is necessary to 
populate the metal-poor regime with more older and more 
massive brown dwarf companions to better constrain this 
region of the parameter space. 
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Table 3. Potential metallicity benchmark binaries 



UCD 


SpT 


Primary 


D 


Sep 


Age 


Mass 


Teff 


loff g 


[Fe/H] 




J - H 


,7 - A's 


Ref 


ID 






pc 


AU 


Gyr 


Mjup 


K 


dex 


dex 


2MASS 


2MASS 


2MASS 




Gl 417 B 


L4.5 


GOV 


21.93±0.21'' 


~2000 


0.08-0.3 


36±15 


1600-1800 




-0.01'^ 


14-57±0.04 


1.10±0.06 


1.88±0.06 


1, 2 


GJ104S B 


LI 


K2Vk 


21.27±0.43" 


250 


<1.0 


55-75 


1900-2200 




-0.07 


13.669±0.160 


0.944±0.258 


1.354±0.140 


3 


LHS102 B 


L8 


M4V 


9.56±0.49'^ 


ISO 


>1.0 


73±10 


1900-2000 


6.0"^ 


0.00 


13.30±0.09 


1.22±0.14t 


1.90±0.07 


4,5 


2M1711 + 4028 B 


L4 + 2 


M4.5V 


21.2±3.9" 


135±25 


1-5 








-1.0^ 


15.00ibO.0G 


0.79*^ 


1,28'^ 


G 


— 1.5 


























HDS9744 B 


LO 


F7IV-V 


39.43±0.4S'' 


2460 


1,5-3.0 


77-80 


2000-2200 




0.17^ 


14.90 


0.90 


1.30 


7, 8 


AB Pic B 


LO-Ll 


K2V 




250-270 


0.03 


13-14 


.„„„+™ 


4.0±0.5 


-O.G4^ 


16.18±0-10 


1.49±0.14 


2.04±0.13 


9 


HD 130948 B 


L0-L4 


G2V 


18.17±0.11'' 


~4S 


<0.S 


<78 


1950±250 




-0.15'^ 


13.9±0.2^ 


0.7±0.2'' 


l.G±0.2'' 


5, 10 


HD 130948 C 


L0-L4 


G2V 


18.17±0.11'' 


~4S 


<0.S 


<68 


1950±250 




-0.15'^ 


14.2±0.2^ 


O.G±0.2'' 


l.G±0.2'' 


5, 10 


fj Cancri B 


L0.5±1.5 


K3III 


91.49±3.35" 


14497-155S1 


2.24-3.50 


G8-72 


1920±100 


5.30-5.50 


0.13±0.06 


17.78±0.06t 


O.G3±0.06t 


1.27±0-06t 


11 


CD-352722B 


L4±l 


MlVe 


21.3±1.4 


G7±4 


0.1±0.05 


31±8 


1800±100 


4.5±0.5 


0.04±0.50 


13.69±0.11t 


0.95±0.16 


1.70±0-13t 


12 


PZ Tel B 


M6 


G9.5 


51.49±2.60'' 


~1G 


0.024±0.003 


62±9 


29S7±100 


4.7S±0.10 


0.05±0.20 


12.26±0.14 


0.39±0.17 


O.S4±0.21 


13, 14 
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5 SUMMARY 

From our age and abundance analysis of the pre-main se- 
quence star PZ Tel A, we have been able to place the low 
mass companion to this star (PZ Tel B) as the youngest 
benchmark object with both an independently measured 
age and metallicity. We find an age for the PZ Tel sys- 
tem older than previous measurements and our new method 
to determine the metallicity of rapidly rotating stars has 
shown that this system is also more metal-rich than pre- 
viously thought. Given the young age of the brown dwarf 
being around 24±3 Myrs, PZ Tel B can provide an anchor 
point to scale evolutionary models at young ages, particu- 
larly when used in conjunction with older benchmarks like 
Wolf 940 B (|Burningham et al.ll2009l ). 

From our age estimates of the host star, we have mea- 
sured bulk properties for PZ Tel B by interpolation within 
a grid of evolutionary models. We find PZ Tel B to be more 
massive and hotter than previous estimates. The mass we 
find (~62±9Mj) places it above the current best estimates 
for the dividing line between planetary-mass companions 
and brown dwarf binary companions. However, we show us- 
ing a basic core accretion derived model that when remnant 
proto-planetary disks contain a high metal fraction, and can 
form large cores (around a few xIO^Mq), core-to-envelope 
ratios of only a few percent can describe well most of the 
current benchmark binaries in metallicity-mass space. 

We ran Monte Carlo simulations to explore the param- 
eter space of our evolutionary models in mass, Toft , and 
log g as a function of metallicity and age. In general we 
found that more metal-rich systems tend towards lower- 
mass companions. Therefore, current models suggest that 
one could increase their sensitivity to lower-mass compan- 
ions in the near infrared by searching around more metal- 
rich primaries. 

We also tested if metallicity has an influence on the near 
infrared colours of brown dwarfs through varying dust opac- 
ity and found a possible trend that suggests that metal-rich 
atmospheres exhibit increased near infrared flux, particu- 
larly in the bluer near infrared bands. We accounted for age 
and mass effects using model corrections and identified a 
possible metallicity trend. However, due to mass/age biases 
in our benchmark sample this apparent metallicity trend 
may be caused by uncertainties in our model corrections, 
and an improved sample is needed. Finally, this paper can 
serve as a blueprint for future analyses of long period ex- 
oplanetary system s that host gas giants like Jupiter (e.g. 
Ijones et al.l I2OI0I ) th at will be directly irn aged by instru- 
ments like the GP I (|Macintosh et al.l l2006f ) and SPHERE 
l|Beuzit et al.ll2006l '). 
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